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Abstract 

We analyze a class of parametrized Random Matrix models, introduced by Rosen- 
zweig and Porter, which is expected to describe the energy level statistics of quantum 
systems whose classical dynamics varies from regular to chaotic as a function of a 
parameter. We compute the generating function for the correlations of energy lev- 
els, in the limit of infinite matrix size. Our computations show that for a certain 
range of values of the parameter, the energy-level statistics is given by that of the 
Wigner-Dyson ensemble. For another range of parameter values, one obtains the 
Poisson statistics of uncorrelated energy levels. However, between these two ranges, 
new statistics emerge, which is neither Poissonnian nor Wigner. The crossover is 
measured by a renormalized coupling constant. The model is exactly solved in the 
sense that, in the limit of infinite matrix size, the energy- level correlation functions 
and their generating function are given in terms of a finite set of integrals. 
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1 Introduction 

Random Matrix Theory (RMT) [[UJ, originally introduced by Wigner, to char- 
acterize the statistical behaviour of the energy levels of nuclei, has found many 
successful applications in various fields of physics in recent years. Originally, it was 
thought that RMT was applicable only to complex systems with many degrees of 
freedom. Hence, it came as a surprise when it was found that it could equally well 
describe simple quantum systems, with very few degrees of freedom, as long as their 
classical dynamics were chaotic. The first evidence of this fact was provided in the 
seminal paper by Bohigas et al p| in which the energy level fluctuations of the quan- 
tum Sinai billiard were analysed and shown to be consistent with the predictions of 
the Gaussian Orthogonal ensemble of RMT. Since this pioneering work, it has been 
checked numerically, on a wide variety of systems, that the local statistical proper- 
ties of a quantum system, whose classical counterpart is chaotic, are well-described 
by RMT. In particular, the nearest neighbour spacing distribution was found to be 
in excellent agreement with the spacing distribution beteen adjacent eigenvalues of 
random matrices. 

In contrast, Berry and Tabor had given strong arguments to justify that, 
for integrable systems with more than one degree of freedom, the nearest neighbour 
spacing distribution of the quantum energy levels should have a Poisson distribution, 
characteristic of uncorrelated levels. This has been confirmed by many numerical 
studies. There now exist excellent reviews on this topic [see || and J7|]. 

However, it is well-known in classical mechanics that purely integrable or purely 
chaotic systems are rare (at least for systems with a few degrees of freedom). For 
most systems, the phase space is partitioned into regular and chaotic regions and 
hence these systems can be referred to as mixed systems. 

An important physical system illustrating these different behaviours, is the hy- 
drogen atom in a magnetic field. The classical system is essentially integrable 
(chaotic) at weak (strong) fields but appears to be mixed at intermediate values 
of the field. This classical behaviour has its counterpart in the energy level statistics 
of the corresponding quantum system, which exhibits a crossover from Poisson to 
Wigner type, when the magnetic field is increased ||. It is, therefore, important 
to find models of random matrices which could describe the statistics of the energy 
levels of such mixed systems. Qualitatively, such a model should be governed by 
a Hamiltonian matrix which is essentially a sum of two parts, one describing the 
chaotic part of phase space and hence belonging to the Wigner-Dyson ensemble of 
the relevant symmetry, and the other corresponding to the regular part of phase 
space. A number of authors have studied models in which block-diagonal GOE 
matrices are weakly coupled by matrix elements also belonging to a GOE ensem- 



ble 12, pj. In this paper, we consider another class of models, first introduced by 



Rosenzweig and Porter [15| to describe the observed deviations from the Wigner- 



and Poisson statistics in the spectra of some transition metal atoms. 
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This class of models is governed by an ensemble of N x N matrices of the form 



H = A+-^G (1.1) 

where A and G are either real symmetric matrices (in the orthogonal case) or com- 
plex self-adjoint ones (in the unitary case). The matrix elements of G are chosen 
to be independent random variables with a Gaussian distribution of unit variance. 
In the context of quantum chaos, G is supposed to correspond to the chaotic region 
of the classical phase space and hence its eigenvalue distribution should obey the 
Wigner-Dyson statistics. In contrast, A should correspond to the classically regu- 
lar region and hence its eigenvalues should exhibit Poisson statistics. It is easily 
seen that the statistics of the energy levels of H depend only on the eigenvalues of 
A. Hence, without loss of generality, the matrix A can be chosen to be diagonal. 
The simplest type of model which can be considered is, therefore, the following one: 
A is a diagonal matrix, whose elements are independent random variables with a 
probability distribution z/(-). Different behaviours can be expected by varying the 
exponent a in (|1.1|). The case a = 1/2, corresponds to a perturbed Wigner-Dyson 
ensemble. It was recently analyzed by Brezin and Hikami They considered the 
case in which G belongs to the Gaussian Unitary ensemble (GUE) and A is a fixed 
diagonal matrix. They showed that the energy level statistics for such a matrix 
ensemble was the same as that of G, i.e., the statistics relevant to the GUE. If 
a > 1, the energy level statistics is expected to be Poissonian. The value a = 1 
corresponds to the crossover regime and for it one expects new statistics. In fact, by 
making a numerical study of this model, Rosenzweig and Porter showed that if one 
chooses the exponent a to be unity, then one obtains energy level statistics which is 
intermediate between Wigner- and Poisson statistics. 

Analytical studies of the model for a — 1, has been done only in the unitary case. 
These studies made use of certain special features of unitary matrices. However, 
the case of the GOE, which one encounters more often, and is technically more 
challenging, had remained virtually unsolved thus far. The only results for this case 
were perturbative ones in the small A limit [|10|j . 

In this article, we develop a technique which can be used to study the spectral 
correlations for the case in which G belongs to the GOE as well as to the one in 
which it belongs to the GUE. We compute the generating function for the average 
value of the product of traces of advanced Green's functions, and the mixed product 
of traces of advanced and retarded Green's functions. All the correlation functions 
of energy levels can be obtained from it, in the limit where N goes to infinity. 

In the case a = 1/2 we show that the statistics for correlations between energy 
levels, on the scale of the mean level spacing, is the same as that of the Wigner- 
Dyson Ensemble. In the unitary case, we hence recover the result of Brezin and 
Hikami M. This matching of our result with theirs is not apriori obvious, because 
the mean level spacing depends on A and we average over the distribution of A, 
whereas they take it to be fixed. 
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In the case a — 1, our result for the generating function (in the infinite iV limit) 
is in the form of a finite set of ordinary integrals. Quite generally, we show that 
the density of states at an energy e is given by v(e), and that all the correlation 
functions are universal functions depending only on the "renormalized" coupling 
constant A = Xu(e). This suggests that in order to make a comparison of the results 
of the model with a given quantum system, it might be plausible to take A to be 
the ratio of p(e) and p reg (e), where p(e) is the classsical Louiville measure of the 
energy surface and p reg (e) is the measure corresponding to the regular part of the 
phase space. 

In order to obtain more concrete results, one needs to evaluate the integrals 
appearing in the expression for the generating functions. This explicit computation, 
which turns out to be a rather lengthy one, has been done in this paper for the 
generating function for the two-point correlation function, in the unitary case. From 
it we can recover the two-point correlation function itself, in the form of certain 
integrals over modified Bessel functions. In the orthogonal case, the generating 
function, even for the two-point correlation function, appears in the form of integrals 
over elliptic functions and we postpone the study of it to a future paper. 

It may be worth giving some hint about the technique used in this paper. We 
basically use integrals over auxiliary Grassmannian variables to compute the average 
over the distribution of the Hamiltonian. However, finally, we evaluate these Grass- 
mannian integrals so as to arrive at a representation in terms of ordinary integrals, 
in the large A" limit. Although similar in spirit to the familiar supersymmetric ap- 
proach, introduced by Efetov M in this kind of problems, our technique is different 
in that we never compute supersymmetric integrals (the only case in which we could 
find supersymmetry useful is when a = 1/2). 

2 Generating Function 

We want to calculate the correlation functions p' n )(ei, . . . , e n ) of the eigenvalues 
Xj of an A" x N self-adjoint matrix H . They are defined as 



pW(e 1 ,...,e n ) = (l[p(e a )), (2.1) 

where 



a=l 



P( e ) = ^E 5 ( e -^)> (2-2) 
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is the local density of eigenvalues at the energy e. The angular brackets will hence- 
forth indicate an average over the probability distribution of H. 

If G s (e) denotes the advanced (s = +1) and retarded (s = —1) Green's function 

G s (e) = ^— -, (2.3) 

e — H — tse 



then 



p(e) = h ™ f— ^tr \G + (e) - G"(e)) \ . (2.4) 
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We will use the following identity to compute [tr G s (e)]/N. 



d det [(e ^ / N)1 N + is (el N - H)} 



de± det [(e+/N)l N + is (el N - H)\ 



isN 



-trG s (e) 



(2.5) 



The symbol Ijy is used to denote the N x N identity matrix. Hence, it is evident 
that the correlation functions of energy levels can be obtained from the generating 
function 

J s 



det [e-(a)/N + is a (e a - H)] ^ 



(2.6) 



*i det [e + (a)/N + is a (e Q — H)\ ' 

where S = {s a }™ =1 , s a G {1,-1} and s±(a) > 0, by taking suitable derivatives of 
it with respect to the variables £-(a) or e + (a). 
In particular, the density of states is given by 



p«(e) 



L 7 RejiAj4 

E- T 0+ I 7T 1 



(2.7) 



= £ +. 



and the two-point correlation function by 



P (2) (ei,e 2 ) 



lim 

T 0^ 



Re 



1 &[Jt- + J} 



2tt 2 9e_(l)9e_(2) 



(2.J 



In this paper, we will consider Hamiltonians H of the form 



H 



12.9) 



where G is an x iV matrix whose elements are independent random variables with 
a Gaussian distribution of unit variance and zero mean. When the matrix elements 
of G are real (complex), the matrix G belongs to the Gaussian orthogonal (unitary) 
ensemble of standard Random Matrix Theory. Since the probability distribution of 
G is independent of the basis, the correlations of the energy levels of H depend only 
on the eigenvalues {aj} of A. Hence, without loss of generality, we can choose A 
to be a diagonal matrix whose elements are independent random variables with a 
probability distribution v(dj). 

2.1 The orthogonal case for finite matrix size 



Consider a mixed system governed by a Hamiltonian matrix of the form ( |2.9| ), 
with the matrix G belonging to the Gaussian Orthogonal Ensemble (GOE). In this 
case, it is convenient to express the generating function J s as the ratio of the square 
roots of the determinants of an antisymmetric matrix and a symmetric one. This 
is because one can cast such a ratio as a product of integrals over real and Grass- 
mannian variables, by making use of the integral identities ( A.l ) and ( |A.2|) . [Note 
that here and henceforth, we suppress the subscript n of the generating function 
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J„.] We can achieve this as follows: We start from the original expression ( |2.6| ) of 
the generating function, which can be rewritten in the form 



with C + and C_ being nN x nN matrices given by 



C 4 



l N -iS®H. 



Here z + , z_ and 5* are n x n matrices with elements 



iV 



S'(aa') = 5 aa ' s a . 
The RHS of ( p.lQ ) can be cast in the form 



(detS)^<^>, 



detM. 



(2.10) 

(2.11) 
(2.12) 

(2.13) 

(2.14) 
(2.15) 

(2.16) 



where S = I2 ® S, and M + and M_ are, respectively, the 2nN x 2nN symmetric 
and antisymmetric matrices given by 



M+ = 1 2 ®C + ■ M_ = T [1 2 ® C_ 
with T being the 2n x 2n antisymmetric matrix: 



(2.17) 



T 







I, 







This leads to the expression 

J s = (detS) N ' 2 { 



"det (TZ_ ®1 N -iT ®HY 1/2 
det (Z+ ®l N -iS®H) 

where the matrices, Z± = I2 ® z±, have the matrix elements 



Z±(pa\p'a') = 8 aa i 8 pp >z±(aa' 



with p = 1, 2. Moreover, 



S(pa\pa!) 
T{pa\p'a) 



o aa i 8ppis OL) 

Sacsn/ipp'), 



(2.18) 



(2.19) 



(2.20) 
(2.21) 
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where 

7 (pp) = ; 7(12) = 1; ; 7 (21) = -1. (2.22) 

We now use the standard trick of expressing the square root of the determinant 
appearing in the numerator of the expression fl2.18| ) as a Gaussian integral over 
Grassmannian variables, and for the one appearing in the denominator, a usual 
Gaussian integral over real variables. [See Appendix A for a summary of some useful 
identities]. In this way the generating function can be written as a superintegral. 

J s = (detS)^ 2 j D<b{eM TZ ^- rsm >)), (2.23) 

where T, Z and S are 4n x 4n matrices defined as follows: 
For a = +, - ; 

Z = Z a 5 aa > ; S = S CT 5 aa > ; T = T a b aa > , 

with 

T + = l 2n , T_ = T, S + = S, S_ = l 2n ; (2.24) 
and $ denotes a supervector 

-(£)■ 

i.e., y? + and are, respectively, vectors with real and Grassmannian elements, 
labeled by 

<p ±i (pa) with i e {1,...,JV}, p G {1,2}, a e {l,...n}. 

At this point, it is important to note a basic difference between the sector of real 
variables corresponding to (p + , and the Grassmannian sector corresponding to If 
in ( |2.13| ) and ( |2.14[ ), all the energies e a are the same and £ a = for all en £ {1 . . . n}, 
then the integrand in ( 2.23|) as well as the measure remain invariant under the 
following change of variables: 

<p +l (pa) = J2 r(pa\p'a') ip' +i (p'a'), (2.25) 

where r is a 2n x 2n matrix satisfying the relation 

t 1 St = S; (2.26) 

and 

<f-i(pa) = U{pa\pa')^_ i {p'a), (2.27) 

where U is a 2n x 2n matrix such that 

U f TU = T. (2.28) 
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The matrices r form a group. If there are p variables s a taking the value +1 and q 
variables taking the value —1, then this group is 0(p, q). Unless q = or p = 0, this 
is a non-compact hyperbolic group. Since the sign of s a corresponds to the type of 
the Green's function (i.e., advanced or retarded), it is evident that this hyperbolic 
symmetry comes into play when we compute the average of a product of p advanced 
and q retarded Green's functions. In contrast, the matrices U belong to the 2n 
symplectic group. We have chosen to extract the factors of s a in the numerator of 
the expression (|2.18 ), of the generating function, in order to emphasize the difference 
between these two sectors. 

Using the expression, ([2.9|), of the mixed Hamiltonian H, we obtain 



J s = (detS)^ 2 [D$ I Du{A) exp \-U,{TZ 



rs<g> A)$ 



x (exp($,(iArS®G)$)), 



where 



N 



(2.29) 



(2.30) 



Dv{A) := J] daj v{ aj ) 
j'=i 

In Q2.29| ) the angular brackets denote an averaging over the Gaussian Orthogonal 
ensemble. 

By integrating over the Gaussian probability distribution of the matrix elements of 
G, we see that 



(exp($, (i^TS ® G)$)) = exp 



where Y is the 2n x 2n matrix defined as follows: 



X 2 



2N 2 



(2.31) 



(2.32) 



with 



Y a (ij) = Y a (jt) 
Hence, 

J s = (det S) N/2 J D<p + J Dip_ J Du(A) exp 
exp \-(<p-, (T_Z_ <S> In- ^- ® A)<p 



,p'=l a=l 



a 



(2.33) 



ip +1 (Z+ <S> In — i$ <8> ^4)¥'^ 



_ ^ 
e 2jv2 



•try 2 



(2.34) 



Since A appears on the RHS of ( |2.31|) only in the form of A 2 , it follows that we can 
always choose A to be positive. We shall make this choice for the rest of the paper. 
We can now decompose tr Y 2 in the following way: 



tr Y 2 = tr (L + S) 2 - tr (L_T_) 2 + 2 (<p_, (T_ ® T)<p. 



(2.35) 
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where the 2n x 2n matrices L a are given by 

L a (pa\p'a') = ¥*i(pa) ¥*i(p'a f ) (2-36) 

i 

and the N x N matrix T is given by 

T ij = J2f+i(P a ) ¥+j(pa)sa- (2.37) 

pa 

We will consider correlations between energy levels around some energy e, on the 
scale of the mean level spacing, i.e., l/(Np(e)). For this purpose we decompose the 
energies e a into 

e a = e + ^ (2.38) 
so that if we define the matrices p + and p_ of elements 

p + {pa\p'a ) = Saa'Spp' [e + (a) + ir a s a ] , (2.39) 
p-{pa\p'a') = 5 aa ,5 pp > [e_(a)s a + ir a ] , (2.40) 

then the generating function can be written as 
J s = (detS^eM*^*-*?-] j Dv(a) j D<p+ j D<p-.e-£** 

.[L'_T_-^] 2 -(v+,(iS®(el N -A)<p+) -(<p_,(iT_® j^.r)p_J 



\2 

x eziv 7 



(2.41) 



2.2 Pseudo— Gaussian transformations 



With the help of Gaussian integrations over auxiliary matrices Q + and the 
exponents in the integrands [on the RHS of (|2.41|) 1 can be reduced to quadratic forms 
in ip + and <£>_. This procedure, which is analogous to the Hubbard Stratonovich 
transformation, yields the following identities: 



x± 



-tr (LiS+^Sp+f = J_ f DQ+e -l tv Ql e -^tr{Q + (L\S+^S P+ )}^ ^ ^ 

d+ J 



and 



with 



-tr(L*_T-i|p_)2 



e2N 

a 



l - J DQ_e-^Q-e->l Q -( Lt - r --$ p -)\ (2.43) 



d± = J DQ ± e-^ tvQ ±. (2.44) 

Our aim is to substitute the above identities on the RHS of Q2.41] ) and change the 
order of integration over ip + , and Q + , Q_. This is because the resulting integrals 
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over if + and turn out to be Gaussian integrals, whose values are obtained by 
making use of the integral identities (|A.l ) and ( |A.2|) . However, the change in the 
order of the integrations imposes a restriction on the form of the matrix Q + . It is 
required to be of the following form [as introduced in [pf] ]: 



Q- 




W5 2 + Q l2 Q 2 i 



Q 



22 



Qn 

tVS 2 + Q2iQi2 



(2.45) 



The block structure refers to the decomposition of the diagonal matrix S [( |2.15| )1 
into p elements equal to +1 and q elements equal to —1, so that Qn is a 2np x 2np 
real symmetric matrix, Q22 a 2nq x 2nq real symmetric matrix, Q12 is a 2nq x 2np 
real matrix and the 2np x 2nq matrix Q 2 \ satisfies 



Q21 



Q 



12; 



(2.46) 



where the superscript t denotes the transpose of the matrix. The variable 5 appear- 
ing on the RHS of ( ^.45| ) denotes an arbitrary positive number. This structure of 
the matrix Q + is dictated by the invariance of the integral over <p + (on the RHS of 
(|2.41|) ) under the pseudo-orthogonal group 0(p,q), where p denotes the number of 
advanced Green's functions and q denotes the number of retarded Green's functions. 



We consider Q_ to be Hermitian : Q[_ = Q_, so for n = 1 it reduces to a real 
number. In addition, we require Q- to satisfy the following relation: 



(Q-T- 



(Q_T_ 



(2.47) 



This constraint is imposed so that the Gaussian integration over the Grassmannian 
variables {(f-j} can be expressed in terms of the determinant of a matrix, as in 
(TOp [see ( [Z3CD - ( 1232]) below]. 

The constraint ( |2.47| ) on the matrix Q_ is imposed so as to ensure the validity 
of the identity (|2.43|) . Indeed the latter identity, involving such a matrix Q_, uses 
the fact that the linear term in Q_, appearing in the exponential, is of the form 
tr(<5— ET_), with E being an antisymmetric matrix. We can therefore write the 
generating function in the form 



(detS) 



N/2 



d + d- 
x exp 



DQ- exp ( --tr \Q H 



.Sp A 



X 



1 

--tr 
2 



K(Q + ,Q-) 



(2.48) 



with 



K(Q+,Q-) 



Dv(a) I Dip + e 



X I_, 



(2.49) 
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where I_ is the integral over the Grassmannian variables {y?-j} and is given by 

I_ = J D V 9„ exp [-(¥>_, Rep-)) , (2.50) 



with 



R := ^(T_Q_) ® Ijv + T_ ® (z(eljv - A) + ^t). 



(2.51) 



The matrix T is defined by Q2.37| ). In the following section we evaluate the integral 
I_ and show that K(Q + , Q_) depends only on the eigenvalues of the matrices Q + 
and Q_. 

2.3 Evaluation of the integral I 

Note that the matrix R, defined by (|2.51| ), is antisymmetric. This follows from 
our choice [ (|2.47|) 1 of (Q_Y_) to be an antisymmetric matrix. Hence, using the 
Gaussian identity, (|A.2f) , we obtain 



I = VdetR. 



(2.52) 



In fact, the constraint (|2.47 ) on the matrix Q_ was imposed so as to obtain the 
above result. 

Moreover, we can show, as follows, that the matrix R depends only on the 
eigenvalues of the matrix Q_: Defining aniVxiV matrix A\\ 



X 2 

A, := i(el N - A) + — T, 



(2.53) 



we write 



detR 



det (T_ 



I at) det 



det (T_ ® 1 N ) det 
det C, 



N 
A 



Q_(g)l N + l 2n ® Ax 
'0 1 N + l 2n ® A x 



(2.54) 



where g_ is the diagonal matrix whose diagonal elements are the eigenvalues of Q_. 



The elements of the matrix C, fl2.54Q , are given by 

A 



C(pa,j\p'a',j') = T_(pp') 



N 



q-{a)8 aa i 8jf + 5 aa >Ai(jj') 



(2.55) 



where q-(a) = q-(pa). This follows from the fact that the eigenvalues of Q_ are 
doubly degenerate (as shown in Appendix B). Hence, the matrix C is antisymmetric 
in the label p and has the form 



C 



D 
-D 
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where D is an nN x nN matrix denned as follows: 



D := ^q- <g> 1 N + l n <g> Ai. 



(2.56) 



In (|2.56| ), g_ denotes an n x n diagonal matrix with diagonal elements g_ Q = g 
This leads to the result 

I_ = detD 



From (g3|) and flT5E| ) it follows that 



n f\ A 2 N 

detD = J] det ljv + i(eljv - A) + — T 

For each a e {1, ... n}, let us define aniVxiV diagonal matrix 

b a := Ijv + *(eljv - -4) 

In terms of this matrix, we can write 



I_ = detD = [] det b a det (l N + j^K 1 ?] 



(2.57) 



(2.58) 



(2.59) 



(2.60) 



For each a, let F a denote a 2n x 2n matrix whose elements are given by 



.'V 



3=1 

It is easy to see that 

tr (b^Ty = tr (SF a y 
for any arbitrary integer j. Hence by ( |A.7| ) we have the identity 



det (l N + ^K 1 ?) = det (l 2n + ^SF a 



(2.61) 



(2.62) 



(2.63) 



Hence from fl2~60l) , ( ggg ) and flgjnf ) it follows that 



N ( \ \ n ( x 2 

I = J] det ^— g_ + i(e - a,)l„,J J] det |l 2n + -j^SF a 



(2.64) 



where denotes a diagonal element of the diagonal matrix A. Using the Grass- 
mannian integral representation ( |A.4| ) for a determinant we can write 



L[det (hn + ^- 2 SF a 



(2.65) 



12 



where tyfi?) is a column (row) vector of length 2n 2 , and 

n In 
a=l p,/3'=l 



(2.66) 



where (3 now refers to the double index (jyy) and B a denotes a 2n x 2n matrix whose 
elements are given by 

Bjfip) = jpS{p)F a {pp), (2.67) 

with <S(/3) = s 7 . The relations ( |2.64|) and ( |2.65| ) yield the following expression for 
the integral I_: 



I 



J e"(*'*) J] det 9- + i(e - a,) I n j 

/ \2 AT 2n n 

k E E v-w03) v+iG90 ^ E [C] , *M *.(^ 

V j=l /3,/3'=l a=l J 



(2.68) 



We can now insert the representation of I_, given by ( |2.68j ), in the expression ( j2.49| ) 
for K, and perform the integration over ip + . This yields 



K(Q + ,Q- 



DVDVe-W D 



>v[a) 



N 



P det ^ q- + i(e - a f )l n 



N 

[]det 

i=i 



5 (i(e - aj)l 2n + ij^Q 



A 2 N 



(2.69) 



where each i?j is a 2n x 2n matrix with elements 

n 

= Ete 1 ),^)*-^ 

a=l J 



n /A \~ — 

E ( + i{el N - A) J tf a (/3) tf a (/3') 



(2.70) 



From the definition (|2.30|) of the measure Du(a), it follows that the expression for 

K(Q + , QJ) involves a product of N identical integrals and hence can be written in 
the form: _ 

K(Q+, Q_) = (det^)-^ 2 / DVDV Z N , (2.71) 



where, writing a for a,j and TZ for 72. j, we define 
A 



Z := J da u(a)det ( — g_ + i(e — a)l n j det 



/ A A 2 N 

Sli(e- a)l 2n + i -fiQ+ + Jp n 



(2.72) 
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That this expression depends only on the eigenvalues q + (pa) of Q + can be easily 
seen by making the following change of variables in the Grassmannian variables 



xM = £<rW>a(/3'') 

P 

where Q is the matrix which diagonalizes Q + , i.e., 

Q+ = Gq+G' 1 , 



(2.73) 
(2.74) 

(2.75) 



and q + is a 2n x 2n diagonal matrix with diagonal elements q + (pa). Hence K takes 
the form 

K = (detS)- N/2 J DxD X e-^ x) P N (2.76) 

where 



P := J da u(a) det g_ + i(e — a)l n J det 
with 



A A 2 ' 
i(e - a)l 2n + + -^ft 



-1/2 

) 

(2.77) 
(2.78) 



We will cast P in a slightly simpler form by making use of the degeneracy of the 
eigenvalues of Q_. Let X± be the In x 2n diagonal matrices of diagonal elements 



and 



x-ipy) = x-(j) := i(e - a) + 4?-(7)- 



We can rewrite P [ ([27770 1 as 



P = da via) 



detX_ 
_detX+ 

The matrix elements of (X^ 1 ^) are 
(X7/ft)(p7bV) = 



1/2 



det | l 2n + ^X + - l K 



1 X a (P7) XaiPl) 



-1/2 



(2.79) 
(2.80) 

(2.81) 



x+(pi) a:- (a) 

It is easy to see that for any arbitrary integer j 

tr (X- l KY = -tr (X^M 



(2.82) 



(2.83) 
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where M is the n x n matrix of elements 



M ma M=it x " (n)x « (i,y) 

p=l 7=1 



(2.84) 



and x_ is a diagonal n x n matrix of elements £-(7). We can therefore write P as 



P 



da via) 



detX_ 
detX, 



nl/2 



det l r 



A 2 M(X 4 

iV 2 



1/2 



(2.85) 



Hence, the multiple integral K(Q+,Q^), defined through ( |2.49| ), is seen to depend 
only on the eigenvalues of the matrices Q + and Q_. It is given by (|2.76| ), with P 
being given by fl2.85|) , above. 

Hence, the generating function is given by 



(det S) 



N/2 



d + d. 



x exp 



DQ + I DQ_ exp ( --tr 



Q+ + i 



X 



-tr 



Q- 



x 



K(Q + ,Q-) 



where 



K = {detS)- N/2 J DxD X e~ ix ' x) P N , 
and P is given by (p.85[) . 

2.4 The unitary case for finite matrix size 

When the Hamiltonian matrix ( [2.9|) , 

A 



(2.86) 



(2.87) 



H 



is such that G belongs to the Gaussian Unitary ensemble (GUE), we arrive at an 
analogous expression for the generating function. 

In this case, we once again proceed from the expression ( |2.1U| ) , where the angular 
brackets now denote an averaging over the GUE. We directly cast this expression into 
integrals over complex and Grassmannian variables by making use of the integral 
identities ( [OP and flO]). This yields 



(ifi + ,C + ip+) -(ip_,C-<fi-) 



(2.88) 



where the matrices C + , C- and S are defined through ( p.ll[ ) - ( |2.15| ). 

The symmetry of the matrices G with respect to unitary transformations allow 
us to use complex variables. There is no necessity to double the dimension of the 
matrices so as to accomodate real variables, as we did in the case of the GOE. 
Hence the label p which appeared in the relations for the GOE, and resulted from 
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this doubling of the dimension, do not appear in the corresponding relations for 
the GOE. All In x In matrices appearing in the case of the GOE are replaced by 
analogous n x n matrices labeled by a single index a. 

Evaluating the average over the GUE in (|2.88 ), we obtain an expression similar 
to ( j2.41| ), with the following changes 



(detS)^ 2 



(detS) N 

In 



(2.89) 
(2.90) 



We notice that in this case the symmetry in the Grassmannian sector is simply the 
unitary one. This implies that the matrix Q_ is now simply an n x n self-adjoint 
matrix. The n x n matrix Q + has the same block structure as before but now Qu, 
with i = 1,2, are self-adjoint matrices and 



Q 



12 



(2.91) 



By proceeding exactly as before, we finally arrive at the following expression for the 
generating function: 



1 



d + cL 
exp 



J DQ + J DQ_ exp (-itr (Q+ + 



ip+S, 
A ' 



|tr(Q. 



P_vi 
A ' 



K{q+,q-) 



where 



K = J DxDxe~ ix ' x) P N ■ 
Here x, \ are Grassmannian vectors of length n, and P is given by 



P :-- 



where 



da via) 



detX_ 
detX 4 



det l n + 



A 2 M(X 4 



E 

7=1 



N 2 X„ 

X*i(l)Xa 2 (l') 



*+(7) 

and X± are diagonal n x n matrices with diagonal elements 

A 



i(e-a) +i— q+(j) 
ar_( 7 ) := i(e - a) + ^-(7) 



(2.92) 
(2.93) 

(2.94) 
(2.95) 

(2.96) 
(2.97) 



respectively. Equivalently, we could also write P in a form analogous to fl2.77| ), 



P := J da v(a) det ( — g_ + %{e — a)T n j det 



A A 2 ' 
i(e - a)l n + i—q + + —11 



(2.98) 
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with the matrix 1Z being defined by ( p.78| ) as before, but the label (3 taking the 
values 1 to n only. 

The above relations, (|2.92 ) - ( |2.94 ), for the GUE, are found to be very similar in 



form to the corresponding relations (|2.48 ) - ( |2.85|) for the GOE. The only difference 



lies in the fact that the expression for P involves square roots of determinants for the 
GOE [see (|2.85|) 1, whereas there is no square root appearing in the corresponding 
relation ( |2.94| ), in the case of the GUE. 



2.5 Supersymmetric formulation for the orthogonal and unitary case 

Up to now we have used Grassmannian variables only as auxiliary variables 
which are to be finally eliminated by integration. We can, however, alternatively, 
cast the expressions for the generating function in a supersymmetric form, where 
Grassmannian and ordinary variables are put on the same footing. If we introduce 
the usual parameter (3 taking the value 1 for the orthogonal case and 2 for the 
unitary one, then the generating function in both cases can be expressed in the 
following elegant form: 



JDQ^e xp (-|Str(Q^ - p^/A) 2 ) F P {Q 
JDQPexp (-iStr(Q/3) 2 



4 = 1 1 T^T. /' 7i, ',J a / n s - ' ' (2-99) 



where 

r ( X 

FpiQ 13 ) = J dav{a) Sdet (e - a)l 4n/p + j^Q 13 J , (2.100) 

the An/ P x An/ (3 supermatrix Q/ 3 being given by 

q9= (qI-2?)< <2 ' ioi) 

and denoting the diagonal supermatrix 




• <l> 



(2.102) 



Here 



S p = 1 2 ®S for 13=1 (2.103) 
= S for (3 = 2. (2.104) 

(2.105) 



The measure of the superintegral is 

DQP := DQ%DQtDQ%_DQL 
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Let us show this in the orthogonal case. We start from the expression ( |2.71| ) for 



K and we note that Z [( |2.72|) ] can be written as 



(det<S) 2 J dav(a)C~ 



where C can be written as 
C = det 

x det 



e - a)l 2n - i-jyQ- 



e-a)l 



2» 



A 

N 



where 



A(H/») 



for a = 1 . . . n, and U is the matrix which diagonalizes Q. 
In this form we see that C is a superdeterminant: 

A 



C = Sdet((e-a)l 4 „ + -Q), 



where 



and 



Q 



N 

Q + Q + - 
Q-+ —iQ- 



Q- 



AU- 1 



Q- 



UA. 



(2.106) 



Q- - At/" 1 ((e - a)l 2n - Q_) Va 



(2.107) 

(2.108) 
(2.109) 



(2.110) 



(2.111) 



(2.112) 



Noting that tr (Q+- Q-+) = ^), we can easily see that the exponential term in 
J 5 , fl2.92p , is indeed of the form 

-istr(Q-f) 2 

The unitary case is treated in the same way. 

We would like to stress that here the supersymmetric formalism gives us only an 
elegant notation. The formalism itself is only useful when we want to mix Grass- 
mannian and real or complex variables. Such a mixing arises in the standard case 
in which the matrix A = 0. This case has been treated by Efetov 0, by using the 
supersymmetric technique. In the large- N limit a saddle point is constructed which 
mixes up the Grassmannian and ordinary variables. The supersymmetric technique 
is also useful in the case where A = 



N and the matrix A is non-zero, where one 
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can show that the energy level correlations agree with the predictions of standard 
Random Matrix theory. However, in the model that we are mainly concerned with in 
this paper, namely the one in which A remains fixed for large N, the supersymmetric 
formulation does not give us any additional advantage. 

2.6 The orthogonal and unitary case in the limit of infinite matrix size 

We will now evaluate the generating function in the limit N — > oo, from which 
the correlation functions can be obtained by taking various derivatives with respect 
to the variables S-{a). It can be easily shown that there is no problem in the 
interchange of the limit N — > oo, with the derivation with respect to £_, and the 
subsequent limit e_ = e + — > 0. 

The result for the generating function both in the orthogonal ((3 = 1) and the 
unitary case (f3 = 2) can be put in the following form: 



J* = 



eX p(^tr (e_-e+) S) 
d + d_ 

x exp ^tr[Q_ - 

PA, 



DQ+Dq_DxDxexp 



A 



(x,x) 



x exp [ \v(e) +Af + A% 



(2.113) 



Here e± are n x n matrices with diagonal elements e±(a), and 

v(a) 



^ f , via) 
p = p( e ) = p / da^^, 
J a — e 



(2.114) 



with the symbol P denoting the principal value of the integral. The matrix Aq is 
given by 

A = i7rtr(q + + iq_)a(q + ), (2.115) 

where cr(q+) is the diagonal matrix whose elements are the signs of the imaginary 
parts of the eigenvalues q + (jry) of the matrix Q + . 



J —oo 



\ det(f + q+) 



3 q+ + iq- 

1 tr 

2 t + q + 



(2.116) 



and 



where 



AS = f_ 



f °° dt f det ( f ~^-) 
oo \det(t + g+) 



det(l -R 13 )? - 1 



J2 j2 x ^"> Xa '^ 



t - iq-(a) p[ ^ t + g+(p7) 



(2.117) 
(2.118) 



19 



and 



Rl a , = — E AaWM 7V ; (2-119) 

This is the main result of this paper. We have expressed the generating function for 
the correlation functions in terms of a finite set of integrals. Hence we have reduced 
the problem of the computation of the generating function, in the limit of infinite 
matrix size, to that of the evaluation of a finite set of integrals. This was our main 
purpose, since, starting from this explicit expression, we can proceed to evaluate 
the physically relevant correlation functions. However, as we shall see, the task of 
evaluating these integrals is non-trivial. Nevertheless, a general conclusion can be 
drawn from this expression by noting that the generating function has the following 
structure: 



J 



exp 



' ip[3 



tr 



E- — e+)S K s ue + , {ur a + 



(2.120) 



where A = Az/(e), can be called the renormalized coupling constant. 

Since the correlation functions can be computed from the generating function by 
using the formula 



1 

2^ 



n 



d 



j cte_(a) 



j(l„ ,-i„) 



=£ + =0 



P (n) (n, 



where 



j(ln-In) _ (TT 



det [s 2 _{a) + (e c 



i det [£%{a) + (e Q - H) 
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(2.121) 



(2.122) 



is positive, it follows that 



-In) _ |_^(In 



-In) I 



and therefore has the structure 



p (n) (ri, . . . , r n ) = v n ft\r x v + a, r 2 u + a, . . . , r n v + a; A), 



(2.123) 



with a = \ 2 f3/2. However, since the correlation functions are translation invariant, 
the RHS of ( |2.1231 ) does not depend on a. 

We shall prove that the density of states p(e) is equal to u(e). We can therefore 
conclude that, on the scale of energy where the mean level spacing is equal to unity, 
the correlation functions are universal, i.e., they depend only on (3 and A. More 
precisely, 



P(e) 



P(e) 



P(e) 



4 n) ( 



,r n ;A 



(2.124) 



Let us now derive eqn. ( |2.113| ). We decompose P, as given by 



da via) 



/ detX_ \ 2 
VdetX + J 



det 



\ 2 M(Xj 



N 2 X_ 



(2.125) 
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into three terms, i.e., P = Pq + Pi + P2, where 



P r X X 

Po := 1 + — I da u(a)ti — — - — - 
2 J X+. 



P\ := da v(a) 



/detX_ 
^detX^ 



P/2 



3 X.-X, 



(2.126) 
(2.127) 



and 



Pi := da via 



Z detX, 
\detX 4 



0/2 



detfl 



A 2 M(X 4 
N 2 X_ 



-1 . 



(2.128) 



Let us first evaluate Pq in the large- iV limit. 



da h>(a)tr 



X_ -X+ 



'n 



da u(a)tx 



(q + + iq_) ((e - a)l 2n + ^q+) 



(2.129) 



since Im g + (7) 7^ 0, the integral on the RHS of (|2.129|) tends to the expression 

— p(e) — i7ip(e)a + (p-y) 

as X — ► 00, where a + (p7) denotes the sign of the imaginary part of the eigenvalue 
9+(P7)- 

Hence, for large N, 



P = 1 + — ^ [ptr (g+ + zg_) + mz/tr (g + + zg_)o"(g + )] 



(2.130) 



where cr(g + ) is the diagonal matrix with elements a + (jyy). In the second term, Pi, 
we make the change of variables e — a = tX/N, so that it reads 



A r+°° / 



det(t — 



/3 q+ + iq- 

1 tr 

2 t + g+ 



. (2.131) 



det(t + q+) J 

The term in the paranthesis is bounded in t and decays like 1/t 2 when t is large, 



since Im q+('j) 7^ 0. Hence, we can use the dominated convergence theorem [OJ] to 
show that if 



sup 



u{t) < 00, 



and u(t) is continuous, then for large N, 



P _ Me) 4/3 



(2.132) 
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The term P2 is treated in exactly the same way as Pi, so that asymptotically, 

\u{e 



Pi. 



N 



Finally, 



lim 



N 



P 



N 



OO 



exp 



'X/3p 



tr (Q+ + iQ_) + \v 



\m +4 + A( 



(2.133) 



(2.134) 



Here and henceforth, we write v for v{e). The expression given by eqn. fl2.113| ) for 
J s is obtained by completing the square in Q + and Q-. 

2.7 The density of states and the average of the product of traces of 
advanced Green's functions 

The only computation which is easy in the general case, is that of the generating 
function for traces of advanced Green's functions. This corresponds to the choice 
s a — 1 for all a 6 {1 . . . n}. 

Let q±(j) denote the eigenvalues of the matrices Q±. For the above-mentioned 
choice of the matrix S, we know that q+(j) = q' + (j) — iS, where q' + (j) is real and S 
is positive. Since q~(j) is also real (and doubly degenerate in the (3 = 1 case) we 
see that the integrands in the expressions for Af and A^ are analytic in the variable 
t in the lower half-plane, and decay like 1/t 2 . We can therefore apply Cauchy's 
theorem to simply conclude that A^ = A^ = 0. In contrast, if we computed these 
quantities for the case of a mixed product of advanced and retarded Green's functions 
(s a = 1, a — 1 . . .p, s a = —l,a — p + l,...n), there would be n — p singularities 
in the lower half-plane, and, therefore, Ay and A^ would be non-zero. 

Hence, it follows easily from ( |2.113| ) that the generating function factorizes as 
follows: 



J ln = exp ( i{p + i-Kv) ^2 [ £ - 



J A- J- 



a=l 



where 



and 



DQ + exp --tr Q + + -*-± ^-(p + tnu) 



J- = — I DQ^ exp 




A 



(p + iiru) 



Hence, we see from the definitions of d± [ fl2.44Q l that J± = 1 and 
J ln = exp [%{p + iitv) ^ [ £ -( a ) ~ £ +( a )] ■ 



a=l 



This implies that 



lim 

iV -> 00 



a=l iV 



) = [—p — mv\ 



(2.135) 



(2.136) 



(2.137) 



(2.138) 



(2.139) 
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which in turn shows that average of a product of the traces of advanced (or retarded) 
Green's functions factorize. In particular, we see that the density of states, p(e), is 
simply given by 

p(e) = u(e). (2.140) 

3 The perturbed Wigner-Dyson ensembles 

We can also look at the situation considered by Brezin and Hikami in the unitary 
case, that is, the situation in which the coupling constant A is of the order of \fN . 
Up to a trivial rescaling, we can simply take A = y/N . We will in this case start 
from the supersymmetric formula ( [2.99 ), which we repeat here for convenience. 

/ZX^expf-iStr^-p/VA) 2 ) , . 

J = " ?— i ~\ 3- 1 

P / J DQ/ 3 exp(-iStr(Q/ 3 ) 2 ) 



where 



A ^ ~ m 



F^) = J dau(a) Sdet ^(e - a)l in/p + -Cfj , (3.2) 
the An/ (3 x An/ (3 supermatrix Q 13 being given by 



and p 13 denoting the diagonal supermatrix 

p" - ( 7 +s * -Z ) • 

The measure of the superintegral is 

DCf := DQ%DQtDQ%_DQt+. 

Following Efetov, 0, let us diagonalize the supermatrix Q, i.e., we define 

Q = VqV-\ (3.5) 

where q is a diagonal supermatrix. Here and henceforth, we suppress the superscript 
(3, unless explicitly required. The measure DQ hence factorizes as follows 

DQ = dfi(u) m(q) dq. (3.6) 

The generating function will therefore be of the form 

Jp = \j dfi(u)m(q)dq exp (--Str^ - ^=^j (Sdet[(e - a)l m + T^]" 1 ^ 

(3.7) 
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where the brackets (■) denotes an average with respect to the measure v(a) on a and 
m = 4n//3. 

Keeping V fixed, we look for a saddle point q, of order yN, of the integral over 
q. We take 

q + = q _ = v/jV (bl m/2 - icSp), (3.8) 

where c is positive and b is real. Since Strg 2 = and Sdet[(e — a)l m + q] = 1, we 
see that the condition for q to be a saddle point is 

" ~(r), (3.9) 



y/N 2 

where 



(e — a)l m + —j= j . (3.10) 



Note that for a function f(a), the symbol (f(a)) denotes the average 

(/(a))= fdau(a)f(a). (3.11) 



The condition (|3.9|) is satisfied if z = b — ic is a solution of the equation (first 
obtained by Pastur jl3l): 

z + ^((e-a + 2)" 1 ) = 0. (3.12) 

Taking into account the contribution of the fluctuations around this saddle-point, 
we see that, if S is the supermatrix 

S = 1 2 ®S, 
then, to order 1 in l/(\/~N), we have 

1 r 

J = - exp(6Str(p)) J dfi{u) m{q) F exp (-icStriVSV^p)) , (3.13) 

where 

F = Jd(5q)ex V (g), (3.14) 



with 



^ = -istr(^) 2 + |str((r^) 2 > + ^ 



2, 



( Str(nfy) ) - (Str(nfy)) 



(3.15) 



However, we see that F and m(q) depend only on b and c, and hence the generating 
function J s has the structure 

J s = exp (i&tr [(e_ - e+)5]) C(6, c) K(ce + , ce_, cr), (3.16) 
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K being the same function as the one in the case where the matrix A is zero, i.e., 
the standard case of Random Matrix Theory. It is easy to see that c = 7rp(e), where 
p(e) is the density of states. Therefore, proceeding as before, we conclude that, on 
the scale of energy in which the mean level spacing is unity, the correlation functions 
are the same as those of Random Matrix Theory. More precisely, 

fe)>fe;ft ;&)-'?W...-.>, w 

p^ being n-point correlation function of the standard Random Matrix Theory en- 
semble, characterised only by (3. 

In the unitary case (j3 = 1), this conclusion agrees with that of Brezin and 
Hikami, [|J, who considered an ensemble with a fixed matrix A. Even though their 
result appears formally to be valid for any set of values of aj (aj being the diagonal 
elements of A), it is clear that at best it is true with probability one with respect to 
some probability distribution on the a/s. In contrast, our result concerns correlation 
functions which are averaged over the a/s. 

In the case a > 1, we take A = l/(iV a_1 ) in the above expressions, and it is 
easily seen that when iV is large, the dominant term in the generating function is 
the generating function for the matrix A. In other words, the statistics is Poissonian 
in this case. 

In order to proceed further, we will now look at the simplest case, that of the 
two-point correlation functions. 

4 Unitary case: The two-point generating function 

There is one major simplification in the unitary case. The integral defining 
A2 [ (|2.117|) 1 can be explicitly evaluated and gives a meromorphic function of the 



eigenvalues {q~j} and {q+j} of the matrices Q_ and Q + , respectively. We will, 
however, only consider the case of the two-point generating function (n = 2), with 
Si = +1 and S2 = —1. The result of the integration over t can be expressed as: 

A 2 = -^[{x 2 + iy 1 )C 22 +{x 2 + iy 2 )C\ x -{x 1 + iy l )Cl 2 -{xi + iy 2 )C n } 
2w z 

TCI 

[detC^ + detC 2 -^ 2 ], (4.1) 



Aw 3 



where 

xj = q +j ; yj = q-j (4.2) 

and 

X\ - x 2 . s 

w = — (4.3) 

We choose Im x\ < and Im x 2 > so that Im w < 0; The matrix C a , with 
a G {1,2}, has for elements 

Cg = X?X* (4.4) 
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while 

Z = C\ X C\ 2 — C\ 2 C 2X + ^11^22 ~~ Cl2^21- (4-5) 

We first integrate over the Grassmannian variables {x'jiX'j} so that if we define 

/ := / DxD X exp(- ]T trC Q + AiM 2 ), (4.6) 

we find that 

/ = F( yi )F(y 2 )g( yi )g(y 2 ) + 3d 2 + d[F( yi ) - g( yi )} [F(y 2 ) - g(y 2 )], (4.7) 

with 

F{y) = l-—( X2 + iy) (4.8) 

g{y) = l + ^ r ( Xl +iy) (4.9) 

(4.10) 

and 

Grouping these results we can express the generating function as: 

j+- = ^/^ + /^-^p[-itr(g + + ^) 2 -itr(g_-^) 2 ; 

x / exp(-Aa[:ri + x 2 + i{yi + 1/2)] + mAz/[-2w + i(y 2 - yi)]) 

xexp (x! +iyi)(x 2 + ^2) (4.12) 

We now integrate over Q_. The integral has the structure 

A = J dQ^F( yi ,y 2 )exp(-hrQ 2 _ + ±tr(p_Q_)), (4.13) 
where F(y 1 ,y 2 ) is some symmetric function of y\ and y 2 . If we define 

Zl := ^ (4.14) 

1/1+1/2 

z 2 := — ~ — , (4.15) 

we can write 

/ z 2 + zi cos if zi sin ipe 1 ^ \ 



Q- 



zi sin ipe ^ z 2 — zi cos y? 
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with ip G [0, 7r] and ip G [0, 27r]. In terms of the variables Z\ and z 2 , we find that 

2ll\ x r+oo /-+00 

A = — 22 a dz i dz 2 z 1 F(z 1 + z 2 ,z 2 - zi) 

x exp(-[^ 2 + z 2 2 \ + ^trp_ + ^tr (p_S)) (4.16) 

and the generating function is given by the expression 

2ttA 



d + d_tr(p„S) ^2A 2 



exp(x^[tr (p + S) 2 -trp 2 _]) 



x 



£>| rfg + exp[-itrg 2 + -^trg + p +1 S]a(^3,w), (4.17) 



where 

* = ^ (4.18) 

and 

^ / \ f r f i t r 2/-. ^TtAz/. 2/1 lH \v 

G a (z 3 ,w) = dz 1 dz 2 z 1 I expL-^H ) - z 2 (l )J 

exp(— niXisw) exp( z 2 Zi — 2Aa[z3 + ^ 2 1). (4.19) 

//• 

This last integral converges if |Im w\ > nXv, a condition that we can impose by 
choosing the free parameter 5 appearing in Q + . 
Now the matrix Q + has the form 



Qa 



Z3 + q — i\/5 2 + h \fh< 



i0 



V Vhe~ i6 z 3 -q + iy/lPTh J 

where 9 G [0,27r], h > 0, q G R and dQ + = 2dz 3 dqdhd9. Since 

hi Q 2 + = w 2 + z 2 with w = ^(q-iVS 2 + h) 2 + h, (4.20) 

we have that 

tr (Q+p+h) = z 3 tr + (g - iv^ 2 + h) tr p+. (4.21) 
Hence, we need to compute an integral of the form 

dq dh exp(-w 2 - itr p + (q - WS 2 + hj) G a (z 3 , w). (4.22) 

- 00 J 

It can be shown that since G a (z 3 ,w) is analytic and bounded in the domain 
Imio< —TrXu, we can change the integral over q into an integral over w along the 



27 



path C := Im w = —b, with b > irXu, by choosing 5 > nXu, so that the integral over 
h can be evaluated. Finally, we obtain 



B 



2i\ 
tr p+ 



f / 11V \ 

J c dww exp(-[w 2 + — trp + ]J (^(23, w). (4.23) 



Hence, 
J+- = 
where 



8n 2 i\ 2 exp(^tr (p 2 + - p 2 _ 
(tr p+)(tr p-.s)d + d- 



p , 11V \ 

^ a I dw w expf — [w 2 +in\h / w+— tr p + ] J K a (w), 



(4.24) 



K a (w) = [ dz 1 dz 2 dz 3 zil expf- V z 2 - ' [(z 3 + iz 2 ) 2 + z\\ - 2\a\z 3 + iz 2 
J v '-r w 



and 
with 



x exp(-*itr (ps) + ytr(p_) - -ytr(p+s)), 
I = A 1 + (z 3 + iz 2 ) 2 A 2 + z\A 3 + [(23 + iz 2 f + z\\ 2 A 4 , 



A, 
A 2 
A 3 
A, 



where 



We now define 



/ iA^ 
(1 + — 
V 2w 


) 4 - 


3A 2 






Aw 6 




A 2 




iA 


"(1 + 


*A 2 


2w A 




2w 3 


2w } 


A 2 




iA 


+ (1 + 


*A 2 


2w A 




2w 3 


2w } 


1 A 4 










16^8' 










A 









M 



f dz x dz 2 dz 3 exp[- VzJ - — \(z 3 + iz 2 ) 2 + z\ 
J 1 ~{ J w 

otz x - y(z 3 + iz 2 ) - u(z 3 - iz 2 ) 



with 



x exp 



u 



2i\ 



tr (P- - P+s) 



y = 2\a - — [tr p_ + tr (p+s) 
t = jtr (p_s), 



(4.25) 
(4.26) 



(4.27) 
(4.28) 



(4.29) 

(4.30) 
(4.31) 
(4.32) 
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so that 



d 



* d 2 A d 2 A ( d 2 9 2 x2 



dy 2 °dt 2 
and the Gaussian integral M is given by 

t 2 



M = 7—r x exp 



A A 



M, 



4(1 + T) +uy - u * 



Since d± = 2tc 2 , we can express the generating function as 



T+- 4t 

.r = exp 

irvt 



d t R, 



where 



R = (wA 1 ) + u 2 (w[A 2 - A 3 ]) + L(wA 3 ) + L 2 (wA 4 ); 
If F(w) is a function of w, we define 

iAu 2 - 



{F(w)) = / dwF(w) exp — w 2 — iw(v + A) 
J c [ 



w 



f +QO , r 2/, 
x / dz exp —z(l H ) - tz 

J-oo 1 w 



and L is the operator 



t 2 a 2 
L = u + W 



We can express i? as a linear combination of the functions 

B n = (w- n ), 

if we note that 

2i 



so that 



LB n = --(n - 2)5 n _x - -5„_ 3 + (u + A)B n _ 2 , 



3=1 



with 



, . 13 2 f A 
ct_i = = 1, a = 2tA, ot\ = A 



iA 3 



A 4 7., vA 3 



a 2 = --iA 6 + — - -wV, a 3 = — - -A + — — + —A + 



2 4 

a 4 = iu 2 A 3 - ^iA 3 , a 5 = -^-A 4 . 

8 4 



4 8 4 16 
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One can also use the recursion formula 



d t B n = -B n - iAd t B n+1 



(4.43) 



to express the generating function in terms of the B n , for n e {—1,4} and d t B 5 . 

One can also show that B n can be represented as an integral over the modified 
Bessel functions. However, we shall not use this representation now, since we do not 
need to compute the generating function itself, but only its derivatives. 

4.1 Unitary case: The two-point correlation function 



In order to compute the correlation function, we start from the expression ( (4.35| ) 
for the generating function and note that when e + = e_, u = 0. Hence, since 
J + ~(e + = e_) = 1, we have that 



dtR 



u=0 



Moreover, 



dJ+~ 



de+{l) 
dJ + - 



de + (2) 



Vt7T 

IT 



—txv + ia, 



-7TU — lOt, 



since these derivatives give the average value of 



1 

"iv tr 



in the limit iV — ► oo. Using these relations, one can show that 



d 2 r 



de + (l)de + (2) 



a 2 + 



v 2 ■ 2A • 

iv H h M 

4 v 



with 



M 



4i 



TXV z 



4(i 



1 d 2 R d 2 R 



+ 



u=0,i;=i= -f (ri —T2) 



•2 9m 2 (9f 2 

Since the two-point correlation function P2( r i; r 2) is given by 

d 2 



2n P2 (r u r 2 ) = Re _ &+(1)ae+(2) 
we see that the unfolded cluster function 



=£+=0 



2 f r i r 2* 



with r = ri — r 2 , 



(4.44) 



(4.45) 



(4.46) 



(4.47) 



(4.48) 



(4.49) 
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is given by 



Y(r) = ^Re {^- V ^-Av-A 2 + mY (4.50) 



where v = nir/A. 

It remains, therefore, to compute M. This is a lengthy computation, which is 
simplified by making use of the recursion formula (|4.43|) and the following relations: 

dv?B n = —iAB n+ i 

d v B n = -iB n _ x . (4.51) 

In this way we get 

M = — £ ^B] + — 2 HdtB,)°. (4.52) 

J=-3 

The (3j are some polynomials of second degree in v, and the superscript indicates 
that the quantities are computed with u = and v — t — mr / A. 

It remains to compute the B n . From eqns. ( 4.37 ) and (|4.39|) , we see that if we 
made the change of variables z = s — iw, when u = and v = t, then we could write 
B n in the following form 

B n = J dww' n J + ds exp (-s 2 - s(v + 2A) - iAs 2 + 2iws) . (4.53) 

When n > 1, we can interchange the two integrals and replace the w integral by a 
contour integral around the origin when s > (because Im w < —A), whereas, if 
s < 0, the w-integral vanishes. In this way one finds that, when n > 1 

B n = 2ni rfs/(s)(2is) n " 1 F n _ 1 (2As 3 ) (4.54) 
Jo 

where f(s) = exp(— s 2 — s(v + 2A)) and 

I n being the modified Bessel function. In order to get this form, we have simply 
expanded exp(— iAs 2 /w) in powers of w' 1 in the w integral. 

When n < 0, we have to do the interchange of the s-integral and the if -integral 
more carefully, and one finds that 

B n = B' n + Bl (4.56) 

where 

B' n = 2ni / dsf(s)(-tAs 2 ) n+1 F n+1 (2As 3 ) (4.57) 
Jo 
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and B'n is some polynomial of degree n in v. Similarly, when n > 1, we have that 



where 



dtB n = iB n _i — R n 

POO 

Rn = n dsf(s)(2is) n F n ^(2As 3 ). 
Jo 



(4.58) 
(4.59) 



We have now at our disposal all the quantities appearing in M, and, therefore, 
from ft4.50|) , the cluster function Y. The final result for Y is given by the following 
expression; 



Y = [°°dsexp (-s 2 - 2As) [Y,a j F j {2As 

+ J™ ds exp (-s 2 - 2As) PjFj(2As 



( ivrs\ 



sin 



A 



where 



and 



A) 
Pi 

ft 
ft 



a = 

Q!x = 
«2 = 

a 3 = 

«4 = 

7A 

47rr 
s 2 A 3 



(« + A) 
4 

s 2 A 
4 

s 4 A - 3As 2 
s s A 2 2A 2 -3 



7r 2 r 2 



8A 2 
7r 2 r 2 



2 -4A 2 



27rr 
5s 4 A 3 



27rr 
4s 6 A 3 



+ s 3 A 2 



irr 

2s 8 A 3 



5A 3 

h 7rr 

27rr 

"A 3 - 7A 


7rrA 


iir 


+ 2 


[1 + 8A 2 


7rr . 


+ T ( 


7rr 





— s 



7rr 



16A 6 



7rr 



(4.60) 



(4.61) 



(4.62) 



This expression appears to be more complicated than the one given in 0, but it can 
probably be transformed into it by using various recursion formula for the functions 



F n , like 



n > 1 xF n+ i = F n _i — nF n and 



dK 
dx 



n+l- 



(4.63) 
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5 Orthogonal case: The two-point generating function 

We start from the equation 



p (2) (ri,r 2 ) 



Re 



d 2 [Jt~ + J 2 H 



2tt 2 \de-{l)de-{2) 



£-=£+=0 



(5.1) 



in order to compute the correlation function p^(ei, e 2 ). Using the fact (established 
in Section p.7|) that 



d 2 j+- 



&_(l)e_(2) 



£-=£ + 

dJ+_ 

aZ(T) 

<9J+~ 
de_(2) 



£_=£-!_ =0 



txv + ip 

7TZ/ — ip 

{-nv + ip) 2 , 



(5.2) 

(5.3) 

(5.4) 

(5.5) 
(5.6) 



we can express p^ 2 - ) (r 1 ,r 2 ) as 
1 



p (2) (r 1; r 2 ) 



2d,d_X 2 7r 2 



z 

DQ+ DQ_ J] [Q_(la|la) + Q_(2a|2a) 



a=l 



where 



with 



x exp(-itr [(Q_ - ^ - */?pA) 2 + (Q + + - ^) 2 ; 



K(q+,q-) = J ,Dxxexp(-(x,x)) exp fA^(e)[~A + A 1 + A 2 ; 



At 



+00 



dt [A(t)\ 



2 2 



1- EE 



?+(P7) +*?-(7) 



and 



A(t) 



2^^! t + g+(P7) 
dt[A(t)]a Idet(l 2 -i2)5 - l] , 

(t-^_(l)) 2 (t-zg_(2)) 2 

n? = in?=i(t+g+(PT)) 
.3.3 



x K{q + ,q_), 
(5.7) 

(5.8) 

(5.9) 
(5.10) 

(5.11) 
(5.12) 



The 2x2 matrix R is given by 



R^ = VrEE^%^- (5-13) 



We recall that q + (jry) and #-(7) are the eigenvalues of Q + and Q-. Two eigenvalues 
of Q + have a positive imaginary part. The other two have a negative imaginary 
part. 

A\ can be expressed in terms of elliptic integrals. The same is true for A 2 , once 
we note that 

det(I 2 - R) = 1 - X with X = tr R - detR, (5.14) 

so that 

8 

det(I 2 -R)* - 1 = ^X n c„, (5.15) 

n=l 

since the X are Grassmannian variables [c n are constants]. One can then integrate 
over the Grassmannian variables, to compute K. We will not reproduce this compu- 
tation here, since we have not analysed the resulting expression for the correlation 
function. It can finally be noted that for higher order correlation functions, hyper- 
elliptic integrals appear. 

Appendix A: Mathematical Preliminaries 

In this appendix we list certain mathematical relations that we make use of later in 
the paper. 

The following Gaussian integrals play a central role in our technique for studying 
mixed systems: 

For a 2n x 2n, antisymmetric matrix M_, we have the relation 

J £V_ e - (v - M -*'- ) = y/detM^, (A.l) 

where f - denotes a vector whose elements ip-j are Grassmannian variables, and 

dip- 



Dtp- :=I1 



-j 



V2 

We use the convention that 



3 -l. 



For a symmetric matrix M + , with a positive definite real part, we have the relation 

J D V+ e-^+' M +^ = [det M + ]" 1/2 , (A.2) 
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where p + denotes a vector whose elements (p + j are real variables and Dp + := 

n^i/v 7 ^) dip +j . 

A similar identity for integrals over pairs of complex variables p + j ) Tp + - ) which is 
valid for any matrix M such that (M + M*) is positive definite, is 

f Vp + e-^+' MLp+) = [det MY 1 , (A.3) 

where _ 

V 2m 

The determinant of any matrix M can be expressed as an integral over Grass- 
mannian variables as follows 

J V(p-e-®-' Mv - ) = det M_, (A.4) 

where Tp_ and p^ denote vectors whose elements, denoted by Tp_- and (p~j, are 
Grassmannian variables; the measure is given by T>p_ = DTp_Dp_, with DTp_ := 
Ylj dtp-j and D(p- := Ylj dip-j. 

The following identities involving determinants and traces of matrices are also 
used frequently: Let A be an (n x n) matrix and B be an (m x m) matrix. Then 

det(A ® B) = (deU) m (det5) n (A.5) 

n m n m 

det(A®l m +l n ® J B) = [] II (a a + b j) = II det i^ 1 ™ + B ) = II det ( A + . 

Q=l jr' = l a=l j=l 

(A.6) 

where a a and bj denote the eigenvalues of the matrices A and B respectively. 
The symbol lj is used to denote the j x j identity matrix. 

• If tr A^ = tr B^ for any arbitrary integer j , then 

det (l n + A) = det (I m + B). (A.7) 

Appendix B: Properties of the matrix Q_ 

The 2n x 2n matrix Q_ is self-adjoint and satisfies the property 

g*_ = -T_g_T_ (B.l) 

where 

T = ( ln \ 

-In 
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We will show that such matrices have doubly degenerate eigenvalues and can be 
diagonalised by unitary matrices, which also belong to the symplectic group. More 
precisely 

Q_ = Uqjfl, (B.2) 

where 

UU ] = U ] U = 1 (B.3) 

and g_ is the diagonal matrix of eigenvalues of Q-. Consider the spectral decompo- 
sition of 

n 

Q- = J2 X 3 P r (B.5) 

3=1 

Then the projectors Pj also have the property ( |B.1| ), which implies that 

Pj(la\l/3) = Pj(2/3\2a). 

hence, 

n 

trP j = 2^P j (la\la), (B.6) 

0=1 

which implies that for almost all Q- (with respect to the Lebesgue measure), Tr Pj = 
2, and, therefore, the eigenvalues of Q_ are doubly degenerate. 
On the other hand, ( |B.1| ) and ( |B.2| ) can be expressed as 

Mq. = q_M, (B.7) 

where 

M = U ] T^{U ] ) t . (B.8) 

However, 

M 1 = -M and M f M = l 2n , (B.9) 



and from ( B.7|) it follows that 

M{pa\qf3) = 8 a>f} [S pA S qt2 - 8 p , 2 S qjl ] (B.10) 

Noting, however, that U can be replaced by Ue lip in ( |B.2| ), where e %v is a diagonal 
matrix, we see from (|B.8|) that we can choose ip a = in ( p.ld|) . In other words 

M= T_ 

which gives the desired property (|B.4|) for U. 
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